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Abstract 

The effect of scattering processes in the continuum on the formation of spectral 
hues in a static atmosphere with an arbitrary distribution of the internal energy 
sources is investigated using Ambartsumian's principle of invariance. Spectral line 
profiles are calculated to illustrate the effect the assumption of the complete redis- 
tribution on atoms and coherent scattering in continuum may have on the emergent 
intensity. The one-dimensional case is considered for simplicity. 



1 Introduction 



In the study of the formation of spectral lines, an interesting problem is the 
transfer of radiation in a medium consisting of centers at which scattering 
occurs for all frequencies. This class of problems includes in particular that of 
the formation of spectral lines in one-dimensional, semi-infinite media. In his 
pionering paper Schuster^ showed that under some circumstances scattering 
of photons in continuum can lead to the formation of emission lines in the 
spectrum. The formation mechanism of emission lines due to scattering in 
continuum is simple. A strong scattering source in the continuum may force the 
photons to be more scattered than absorbed in the certain part of the spectrum 
and will therefore lead to a decrease in the source function in the continuum. 
Suppose that the radiation field in a spectral line is in local thermodynamical 
equilibrium (LTE) with the local medium so that the line photons have much 
more chance of being absorbed and destroyed in the medium than photons in 
the continuum. The increase in the efficiency of scattering causes the radiation 
field in the continuum to decrease in the meantime, thereby keeping the flux 
in the spectral line unaffected because of the thermal radiation of the medium. 
Thus the spectral line will appear in emission. 
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2 Spectral lines formed in an isothermal medium 

Recently Israelian and Nikoghossian^ investigated the diffusion of isotropic and 
coherent radiation in a semi- infinite medium by the probabihstic approach. 
The probabilistic method is based on the concept of a probabihty of quantum 
exit from a medium and has been successfully applied to treat many classical 
problems in a theory of radiative transfer^ . The principle of optical reversibil- 
ity, together with the principle of invariance, allows very complex problems to 
be handled when the probabilistic approach is employed. 

Let us consider a one-dimensional, semi-infinite, isothermal medium with a 
distribution of energy sources: 



where u{x) = (1 — X)a{x) + j3. Here we denote a; as a non-dimensional fre- 
quency {x = 5v/5vd: where Svd is a Doppler width), a{x) as the normalized 
absorption coefficient profile in the line, i?[T(r)] is the Planck function, which 
depends on the optical depth, r, relative to the central wavelength of the fine 
through the temperature, T, A is the probability of photon re-emission, (3 is 
the ratio of the absorption coefficient in the continuum to that at the cen- 
ter of the spectral line. We have also introduced^ a parameter 7, the ratio 
of the scattering coefficient in the continuum to the absorption coefficient at 
the line center. Let us assume that the medium has an infinitely large optical 
thickness so that the addition of a layer of small optical thickness will not 
change the reflection coefficient of the medium. This principle of invariance 
of Ambartsumian^'^ allows one to obtain the following equation: 



where v{x) — a{x) -|- /3 -|- 7, p{x,x') and r{x,x') are the reflection and re- 
distribution functions, respectively. The reflection function, p{x,x'), has the 
following probabilistic meaning: if a photon of frequency x' is incident on the 
medium, then p{x,x')dx is the probability that, after multiple scattering, a 
photon is reflected from the medium with a frequency in the interval x,x + dx. 

Hereafter, all integrations over x will be carried out from —00 to +00. In 
general, there can be several sources of absorption/ scattering in the medium, in 
which case r{x, x') must be considered as a sum of redistribution functions due 
to different centers (like atoms, electrons, dust particles, molecules, plasmons. 



e{T,x)^u{x)B[T{T)l 
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etc.). Let us also add pure scattering centers in the medium and denote their 
redistribution function as rg{x,x'). If these centers are free electrons and the 
scattering is incoherent^ then 

oo 

re{x, x') = 5i/i3-V-°-^(e-^' - 2xTT J exp{-u^)du). (3) 



In case of the coherent scattering we have: 

reix,x')^j5{x-y). (4) 



The total redistribution function in the case of the complete redistribution of 
photons on the atoms and the coherent scattering on the electrons will be 

r{x,x') — a{x)a{x') + —6{x — y). (5) 

A 



The problem of determining the radiation field in a medium with or without 
energy sources under broad assumptions concerning the elementary scatter- 
ing process can be simplified if the redistribution function for atoms allows 
presentation in the form of a bilinear expansion^'^'^: 

CO 

r{x,x') = ^Akak{x)ak{x'), (6) 

k=0 



where 



1 1 _ 2 

ak{x) = ^ H2k{x) 



2k + 1 7rO-252fcJ(2A0! 



and H2k{x) are the Hermit polynomials. Note that with only the first term 
taken in the expansion (6), we arrive at complete frequency redistribution, i.e., 
there is no correlation between the frequency of the incoming and the scattered 
quantum in the observer's frame of refernce. Such a representation can be 
obtained by expanding the function r{x, x') with respect to its eigenfunctions 
over (—00,00). In general, the function r{x,x') can be replaced by an finite 
sum of its expansion with respect to its eigenfunctions over (—00,00): 



N 
k=0 



r(x,x ) rjv(x,x ) = 2^ , (7) 
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where the functions Xfe(^) ^^e the normahzed solutions of the equations 

oo 

Xk{x)^Xk j r{x,x')xk{x')dx'. (8) 



For fixed N, this partial sum gives the best mean-square solution approxima- 
tion in the (x', a;)-plane to the function r{x',x) among all possible represen- 
tations. In the case of representation (7), the accuracy of the fulfillment of 
condition (8) agrees with the accuracy of the approximation of the function 
r{x,x') by the partial sum risf{x,x'). Different methods (like the method of 
least squares, the method of moments, etc.) exist for the approximate con- 
struction of the eigenfunctions of the kernel r{x, x') and we shall not dwell on 
this point. 

Suppose that the function re(x, x') can be expanded as well: 



re{x,x') = ^ BmUJm{x)i^m{x')- 
m=0 



(9) 



Then in a more general case we will have a sum of functions (6) and (9). If 
the sum of representations (6) and (7) for the redistribution function, r{x,x'), 
is used, then the integral equation (2) can be rewritten as: 



p{x,x') = ^ 



Lk=0 



'v{x) + v{x') 



i=Q 



'v{x) + v{x') 



(10) 



where the functions (pk and ipi{x) satisfy the set of functional equations provid- 
ing a generalization of Ambartsumian's functional equation for non-coherent 
scattering: 



(Pk{x) = ak{x) + J ak{z)p{x,z)dz, (11) 

— oo 
oo 

'ilji{x) ^ u;i{x) + J uji{x)p{x,z)dz. (12) 



The solution of this system provides a reflection function, p{x. x'), which allows 
the evaluation of the profiles of spectral lines (since R{x) = J p{x,x')dx'). 
Determination of the profiles of absorption lines in an isothermal atmosphere 
is equivalent to that for diffuse reflection of radiation from a semi-inflnite 
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medium^. The reflection and absorption profiles in this case are related by 

R{x)+R{x) = l. (13) 



This relation enables one to reduce line-profile calculations for an isothermal 
atmosphere to the evaluation of the reflection function, p{x,x'). 

In this article we present numerical results for the simplest case of complete re- 
distribution of photons on atoms and a coherent scattering on free electrons. If 
last term in the equation (5) describes a coherent scattering in the continuum, 
then it can be written in the form 

^5{x - x') ^^Yl'^k{x')ak{x), (14) 
^ ^ k=0 



using the set of the ak{x) functions given in the expansion (6). 



3 Spectral line formation in a non-isothermal medium 

In the case of a non-isothermal medium we shall assume that the source func- 
tion can be presented in the following form: 

B[T{t)] = E -^((^^r^ (15) 

n=0 ^■ 

where Bn is the n^^ derivative of the Planck function. Let us denote by 
Y{t,x',x) the probability that a photon with a frequency x' moving in the 
medium at an optical depth r will leave it with a frequency in the range 
{x,x + dx). Similarly, let P{t,x',x) be the same probability but for photons 
absorbed at an optical depth r. Application of the invariance principle and 
simple physical reasoning results in the following equations: 



dY{T, x', x) 



-v{x)Y{t,x',x) + J P{0,z,x)a{z)Y{T,x',z)dz (16) 



dr 

A 



— oo 

oo 



A A f 

a{x)P{Q,x',x) = —r{x',x) -\ — / r{x',z) p{z,x) dz (17) 
2 2 t/ 



These equations can be easily derived when (3=0 (Sobolev^) and (3 ^0 (Nikoghos- 
sian and Haruthyunian^) . 
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The intensity of the emergent radiation, 7(0, x), with the source function given 
by (15) can be presented in the form 



I{Q,x) = J2BnW,x), (18) 



n=0 



where 

oo oo 

4(0, x)= I u{x')dx' I fn(T)Y^r,x',x)dr (19) 

-oo 

and /„(t) = (/3T)Vn!. 

The spectral hne contours for the mentioned distribution of the energy sources 
can be calculated with the formulae^ 

oo 

v{x)Rn{x) = I a{x')P{0,x',x)Rr,{x')dx' + f3Rn-i{x). (20) 



Note that 7 appears only on the left-hand side of this equation. The final 
profiles are 

Rix) = ^^^^^iM. (21) 



Note that the function p{x,x') in (17) is the reflection function from a semi- 
infinite isothermic medium. Since we have presented this later in the form of 
eq. (9), we can write, instead of eq. (17), 

2 00 00 

-a{x')P{Q, x\ x)^Yl Akak{x')ipk{x) + ^ B^ujm{x')'il)m{x)- (22) 

^ fc=0 m=0 



Let us multiply eq. (20) by ak{x) and integrate with respect to x over the 
range (—00 cxd). Repeating the same operation for uim{x) yields, for — 
J a{x')Rn{x')dx, the following system of the algebraic equations: 



Rn=i: AkR% + E BmR'^O,,^ + 4"-'^ (23) 
fc=0 m=0 

00 00 

i?- = E AkRl^J.m,k + E BmRnOm + Ci""'^ , (24) 
fc=0 m=0 
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where we have introduced the notation 



J V{X) J v{x) 



-oo —oo 

CO , , CXD 



-oo — oo 

oo , , , , , oo 



J vix) J V[X) 



Apparently the question is reduced to solving the system of algebraic equations 
(23-24). It is obvious that we can add different absorption/scattering centers 
(atoms, molecules, etc.) to the medium, and if their redistribution functions 
allow bilinear expansions, the problem can be reduced to a system of algebraic 
equations. For a medium containing two kinds of absorbing/scattering centers, 
one needs to construct and tabulate functions (pki^) and iprn{x)- Once this has 
been done, one needs to compute the constants ^k, o'k,m, S^^K l^m,ki Omi and C,!^\ 
after which the line profiles are evaluated from eqs. (20) and (21) using recur- 
sion relations. This discussion clearly remains valid for more than two kinds 
of absorbing/scattering centers in the medium (provided that their redistribu- 
tion functions allow a bilinear expansion) and for more general assumptions 
relative to the geometry of the medium. 



4 Discussion 



It is known that coherent scattering in a medium gives the source function 
a non-LTE character. This fact alone makes it possible for spectral lines to 
appear in emission, even if the line source function does not deviate from 
LTE^. The mechanism requires a low temperature gradient, a strong spectral 
line and an LTE process of line formation. The criteria derived for the appear- 
ance of the lines in emission in a non-isothermal stellar atmosphere are much 
more stringent and it is possible that this mechanism cannot work in a real 
stellar atmosphere. Another difficulty in producing emission lines in a normal 
atmosphere is the small cross-section of the Thomson scattering. However, in 
the general case, Thomson scattering on free electrons is not necessarily the 
only nor the most effective mechanism among the scattering processes which 
may exist in astrophysical plasmas. 

We have computed line profiles from a semi-infinite, one- dimensional atmo- 
sphere assuming complete redistribution on atoms and coherent scattering on 
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free electrons. For the isothermal and non-isothcrmal (^^=1) cases we have 
numerically solved the systems of equations (10-12) and (23-24), respectively. 
Problems connected with the method of approximation for the solution of the 
system of equations of the form (10-12) are discussed in detail in the work of 
Haruthyunian and Nikoghossian^. 

Figures 1 and 2 illustrate the absorption/emission profiles for different values 
of 7, A, and f3. The profiles for the isothermal and non-isothermal {n=l) 
media are shown in Figs. 1 and 2, respectively. To compute the S„ coefficients 
in the sum (15) we have considered a hypothetical fine at 6000 A, a stellar 
atmosphere with Tcs = 40000 K and a gray-atmosphere temperature law. We 
can see that even in the case of a non-LTE source function in the line (i.e., 
A <1) one can have emission lines depending on the efficiency of scattering and 
absorption in continuum. The five parameters involved in this problem are the 
scattering (7) and the absorption (/3) in continuum, the scattering (A) and the 
absorption {<y{x)) in the line and the temperature gradient in the medium. The 
addition of scattering to the atmosphere can raise or lower the local continuum 
depending on the location of the spectral line, the efficiency of scattering and 
a temperature gradient Therefore, the variations of continuum in Figs. 1 
and 2 have only a local character. We have already seen^° that the variation 
of the continuum due to the scattering is nothing else than a redistribution 
of emergent radiation demanded by the condition of radiative equilibrium and 
the conservation of integral fiux. Our Fig. 2 clearly demonstrates that the 
conditions for emission are much more stringent in the non-isothermal than in 
the isothermal atmosphere. It is also possible to have a situation in which the 
wings (formed in LTE deep in the atmosphere) of a strong line will appear in 
emission while the line core is in absorption because of the strong non-LTE 
effects. Interesting effects may appear when one studies a wing and a core of 
the spectral line separately, considering depth and frequency dependence of 7, 
(3, and A. 
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Fig. 1. Spectral line profiles in the case of complete frequency redistribution on atoms 
and coherent scattering in continuum for an isothermal atmosphere. Computations 
have been done for 7=0.001 (solid line), 7=0.1 (dash-dotted line) and 7=1.0 (dotted 
line). Values of /? and A are indicated. 
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Fig. 2. The same as Fig 1. but for a non-isothermal atmosphere with n = 1 (see 
text). 
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